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ABSTRACT 

Turbulence is a key element of the dynamics of astrophysical fluids, including those of interstel- 
lar medium, clusters of galaxies and circumstellar regions. Turbulent motions induce Doppler shifts 
of observable emission and absorption lines and this motivates studies of turbulence using precision 
spectroscopy. We provide high resolution numerical testing of the two promising techniques, namely, 
Velocity Channel Analysis and Velocity Coordinate Spectrum. We obtain an expression for the shot 
noise that the discretization of the numerical data entails and successfully test it. We show that nu- 
merical resolution required for recovering the underlying turbulent spectrum from observations depend 
on the spectral index of velocity fluctuations. Thus the low resolution testing may be misleading. 
Subject headings: 



1. INTRODUCTION 

As a rule astrophysical fluids are turbulent and the 
turbulence is magnetized. This ubiquitous turbu- 
lence determ ines the transport properties of interstellar 
medium (seelElmegreen fc Falgaron i[i99l,|Stutzkl|200i|, 
iBalesteros-Pered es et a ljl2006l) and intracluster medium 
(see ISunvaev. N orman fc BrvanI l2003l lEnsshn fc Vog^ 
120061 lLazarianH2006ll . manv properties of So lar and stel- 
lar winds (see iHartman fc McGrcgol [ToMl etc. One 
may say that to understand heat conduction, propa- 
gation of cosmic rays and electromagnetic radiation in 
different astrophysical environments it is absolutely es- 
sential to understand the properties of underlying mag- 
netized turbule nce. The fasc i nating processes of star 
formation (see iMcKee fc Tal] [20021 lElmegreenI [20021 
Mac Low fc Klessenll2004[ ) and interstellar chemistry (see 



2003 ) and also with absorption effects in iPadoan et al.l 



Falgarone et al.l (|2006( ) and references therein) are also 



intimately related to properties of magnetized compress- 
ible turbulence (see reviews by Elmegrccn fc Scalo 2004]) . 

We should stress that while the density fluctuations 
are an indirect way of testing turbulence, the most valu- 
able information is given by the velocity field statistics, 
encoded in spectrometric observations. To recover it we 
need to account for mapping of a velocity field in real 
space to position-position-velocity (PPV) observational 
data cubes. 

We have several options here. For subsonic or 
slow supersonic turbulence we can use techniques 
of velocity centroids, wh ich allows us recove ring of 
velocity spectrum (Esqu ivel fc LazarianI l2005f ). As 
it was shown there, this technique does not work 
for turbulence with high Mach number. For such 
turbulence Velocity Channel Analysis (VCA) (see 
[Lazarian fc Pogosyan 2000), which deals with 2-d slices 
of PPV, and Velocity Coordinate Spectrum (VCS) (see 
Lazarian fc Pogosva nll200dl Lazarian fc Pogosva a [20061 . 



Chepurnov fc Lazar ian 200(|) should be used. 



VCA has been already used with both atomic hydro- 
gen and molecular line data to r ecover the und erlying 
turbulence spectra (see review bv lLazarian[ 12006( 1. The 
technique ha s been tested with nu meric s for optically 
thin case in ILazarian et al.l (|2001h and lEsauivel et al.l 



20061) . However, rather strong shot-noise was observed 
in the synthetic channel m aps. For instance, ri umeri- 
cal simulations for VCA in lEsauivel et al.l (|2003f ) show 
that discretization over line of sight may distort substan- 
tially the spectrum at high wavenumbers. It was also 
shown, that its impact can be minimized, when veloc- 
ity slice is made thickest possible for the thin (velocity- 
dominated) regime. The shot-noise was also observed 
to decrease with the increase of the cube size used for 
analysis. Nevertheless, the exact origin of it was some- 
what unclear. This effect seemed to confuse other re- 
searchers who started to wonder about the practical util- 
ity of the technique with th e real spectroscopic data 
([Miville-Deschnes et al.ll2003D . 

The VCS is a more recent devel opment. Although the 
formu lae for it were obtained in ILazarian fc Pogosvanl 
20001). the actual app hcation of it has just started 
Chepurnov et al.[ l2007l ). No extensive testing of the 



technique has been performed before this work, as far 
as we know. 

In this paper we deal with the numerical testing for 
both VCA and VCS techniques. In particular, we subject 
the origin of shot noise to scrutiny in this paper. 

In what follows we provide some facts about the VCA 
and VCS in Sect. [2l analyze shot-noise in Sect. [3l de- 
scribe the results of numerical testing in Sect. [4| and 
provide the summary in Sect. [6l 

2. BASIC FACTS 

We show in Fig. [1] the essence of the practical data 
handling with the VCA. The data in the PPV data cube 
is being analyzed in the slices of AT^ thickness. The 
images of the eddies can have an extend both larger 
and smaller than AF . The former case is termed by 
ILazarian fc Pogosvanl (|2000[ ) "thin slice" , the latter is 
termed "thick slice" . 

In Fig. [21 we illustrate the essence of the VCS tech- 
nique. Depending on the resolution one may or not re- 
solve the spatial extent of the eddies under study. 

Below we briefly o verview the main a n alytica l 
results, obtained in ILazarian fc Pogosvanl (|2000| ). 
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iLazarian fc Pogosvanl (l2004f). iLazarian fc Pogosvanl 



( 20061 ) and IChepurnov fc LazarianI ()2006l ). which are 
relevant for this paper. 

Our calculations are based on the following expression 
for a spectral line signal, measured at velocity vq and 
given beam position: 



wir) dfe{f)fivr{7^ + «mf) - «o), (1) 



where e is normalized emissivity, / is convolution be- 
tween channel sensitivity function and Maxwellian dis- 
tribution, defined by temperature of emitting medium, 
and w is a window function defined as follows: 

w{f) = ^Wb{c^,e)weir), (2) 

where wi, is an instrument beam, is a window function 
defining the extend of the observed object. 

If we take a Fourier transform of S and correlate this 
quantity taken in two directions 1 and 2, we have the 
following measure, which can be used as a starting point 
for t he mathematical formulation o f both VGA and VCS 
(see IChepurnov fc LazarianI (|2006f) for details) : 



Ki2{K) = (Si{ky)S2 (K) 



fc2 

exp ( --^D^ir) - ikybz 



where 



Wi2{f) = / wi{T')w2(r' + r) dr' 



(3) 



(4) 



and is a velocity structure tensor projection: 

D,{r^-f') = ({v,if)-v,{r')f), (5) 



Expression [3] is true provided that the following con- 
ditions are satisfied: 1. fields involved are homogeneous, 
2. velocity is Gaussian, 3. absorption effects are neg- 
ligible^, 4. beam separation as well as beam width are 
small enough to neglect the difference between Vr and Vz 
(we consider z-axis to be a bisector of the angle between 
beams). We also assume that 5. ^^'^^(r) depends only on 
z and admits a linear approximation with angular coefh- 
cient b. The relative importance of these conditions and 
the way to avoid some of them wi thin the VGA and the 
VGS frame work are discussed in ILazarian fc Pogosvanl 
dlOOi) and 'Lazari an fc Pogosvanl (|2006'). However, for 
the sake of simplicity, we do not discuss the most general 
case here. 

If we set wi2 = Wb,a(,R— Ri2)w^,a{z) (i?i2 is a beam 
separation in picture plane, Wh,a and w^a are auto- 
convolutions of Wb and w^), which corresponds to the 
case of remote emitting object, and take Fourier trans- 
form over J? i 2 , we have a power sp ectrum, designated in 
VGA as Ps ILazarian fc F'ogosvan|[2000l ): 



Ps{K,K) = f{K)Wb\K) I w,{z)dr 



e-'(™'=^^)C,(r)exp(-^i^,(f) 



(6) 



^ In this case analysis is more simple if we use the Fourier space. 
Absorption effects can be described only in real space, which im- 
plies u s ing of structure functions ins tead, see lLazarian fc Pogosvanl 
II2004I ) . ILazarian" fc PoeosvanI ll2006h 



It is convenient to introduce P3 as a 3D power spec- 
trum of a PPV data cube with assumption of infinite res- 
olution over V and R, which, for instance, corresponds to 
the conditions for the unprocessed synthetic PPV cubes: 



(7) 



P:i{K,K) = {^\S{K,K)\''j= J We{z)dr 

•e-'(™'=-)C,(rOexp(^-|i^,(f) 

In this case P^J?, = p{ky)m^{K)P3{K, fc„). 
The VGA predicts the behavior of the quantity 



P2iK)= / dkyPsiK,K) 



Wb\K) / dk,f{k,)P3{K,k,), 



(8) 



which can be easily determined from observations. For 
the channel which effective width^ is small enough (see 
Fig. |T]) P2 is velocity dominated regime with the slope 
slope {9 ~ av)/2, provided that the density has a steep 
spectrum, i.e. > 3. the later can be estabhshed 
from column density measured with infinitely wide slices. 
Then the spectum is density dominated, and the slope 
is a^. Here and are the slopes of velocity and 
density power spectra . For a more detailed analysis see 
ILazarian fc Pogosvanl (|2000( ). while the predictions for 
different variants of velocity-dominated regime are sum- 
marized in Tab. 

The one-dimensional spectrum Pi is the subject of the 
VGS studies. It corresponds to the case when the two 
beams, involved in K12, coincide, i.e. 

Pi(fc,) = ifi2(fc.)|^^^^o = P{ky) I wnir^dr 

■C,ir)cxpi-^Dz{r}^ik,bz 

In terms of P3 the one-dimensional spectrum can be 
written as follows: 

Pi(fc,) = P(fc,) / dKm\K)Pz{K,ky)- (10) 



As it was shown in ILazarian fc PogosvaiH (|2006[ ) and 
urnov fc Lazaria nl (|2006D . here we also have two 
regimes of spectum. They depend on beamwidth, see 
Fig. O for illustration. For high resolution mode {k^ is 
less than velocity variance on the beam scale) the slope 
of Pi is 2/(a„ — 3), otherwise, in low resolution mode it 
is 6/{av — 3) (steep density spectrum is assumed for the 
both cases) . More complete list of the predictions for Pi 
is presented in Tab. O 

If we compare Eqs. [8] and [10], we see that the VGA 
and the VGS provide complimentary approaches and are 
related to the same quantity, P3. Their properties must 
be closely related therefore. 

^ It is approximately equal to \/(^^ch (■^'"t)^), where Sv^i^ is 
channel width and vx is thermal velocity of emitting atoms. 
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3. SHOT-NOISE 



In lEsQuivel et alj (j2003[ ) it was shown, that the quality 
of the VGA simulations strongly depends on the number 
of points along the line of sight. The underlying interfer- 
ence was termed shot-noise to reflect is stochastic nature: 
insufficient statistics for the signal in spectrometer chan- 
nels. We revealed the same type of noise for VCS. In this 
case, Pi saturates to a constant level, which decreases an 
available interval of power law (see Fig. O right). Below 
we address the condition needed to avoid it. 

To have a realistic spectral line we need a velocity field 
to be resolved over line of sight well enough (see also 
llnogamov fc SunvaevI [20031 ) . This means, that the typ- 
ical scale over z, for which velocity stays within a sin- 
gle channel near its extremum, should have enough dis- 
cretization points, because such extrema give the domi- 
nant contribution to the spectrometer signal. 

Below we estimate the required resolution. To do so, 
we need to bind the correspondent typical scale over line 
of sight with the channel width. The natural choice is 
calculating velocity dispersion (which itself can be associ- 
ated with the channel width) for the wavelengths shorter 
than this scale. As mentioned above, the contribution of 
the velocity profile shapes other than extrema, provided 
by the remaining harmonics, can be neglected. There- 
fore, the scale found by reverting of this expression is 
what we need. 

In other words, we divide the 1-d velocity spectrum 
domain into the intervals separated by the scale in ques- 
tion La, and, as spoken, assume that the harmonics in 
its low-wavenumber part provide an extremum of veloc- 
ity. For the higher harmonics, in order to place the bulk 
of the signal inside one channel, we should demand, that 
the correspondent variance (j{L/^) is about half a channel 
width: 

.(..)»^. ^ ^ ^ (U) 

where La is a scale, at which velocity stays within a 
single channel near its extremum, Na is a number of 
channels per velocity variance, and cr(L) is a velocity 
variance for scales smaller than L, expressed through 1- 
d power spectrum'^ Pi^: 



a2(L) = 2 J Pu{k)dk. 



(12) 



If a 3-d power spectrum has a cutoff at fco = 2tt/L, 
where L is an injection scale, Pid has the following form 
for the isotropic 3-d spectrum of 

1 r fc-(""-2), \k\ > ko 



Pi 



Id 



2 I fc-(""-'\ |A:| <A:o 



Then, having solved Eq. [TT] for La we get 
2N„ 



La 



L . 



(13) 



(14) 



If TiA is a number of points per La, then the total 
number of points along z is as follows: 



riA- 



2N„ 



L y^Ja^ — 2 



(15) 



Should not be confused with Pi , that refers to the PPV space 



where L^ is the size of emitting structure, riA is deter- 
mined from simulations and is experimentally found to 
be around 6. 

The requirements for Nz are rather demanding. If we 
have taken our typical settings: 128 channels per triple 
line width and Ls/L = 2, then for = 4: Nz there 
should be at least 1.1 • lO"' points, for = 11/3 the 
corresponding number is 4.3 • 10^, while for a„ — 3.5 the 
number is as big as 1.7 • lO''. Therefore it is not surpris- 
ing that without this knowledge so me researchers could 
get p uzzled by their results (e.g. Miville-Deschnes et al.l 

mm . 

Condition |15| guarantees, that the spectrum is clear 
from shot-noise for all available fci,'s, i.e. up to 2'k/IS.v, 
where Aw is the channel width. This estimation refers to 
the case without smoothing along angular coordinates, 
i.e. for the high-resolution regime. Low-resolution VCS 
spectrum drops faster and therefore the requirements are 
even more restrictive. 

4. RESULTS 

VCS spectra, illustrating condition [15] are presented 
in Fig. [3l The velocity spectral index decreases from 
the top to the bottom (the values are 4, 3.78 and 3.56). 
According to our estimates above, this requires a drastic 
increase of number of points. Indeed, Nz is, respectively, 
2^**, 2^^ and 2^^ . In the left column with the shot-noise 
clearly present Nz is 4 times undervalued, in the right 
one it is increased by factor two. 

Examples of restoration of the velocity spectral indices 
from simulations are presented in Tab. [T] They show 
good agreement, with the exception of somewhat out- 
falling value for a„ = 4. Values of Nz have been chosen 
regarding Eq. [15] . 

Knowing the requirements on Nz , we generate the VCS 
and VCA spectra (see Fig. |4l with = 4). But for 
shallower velocity spectrum satisfying such conditions for 
VCA may be difficult, even if we use 2-d data instead of 
3-d data cube. For example, the condition [15] for the 
VCA simulations at = 11/3 implies a square array 
with the side 420000, which is difficult to implement at 
the moment. If we decrease this size, the VCA spectrum 
gets completely distorted, see Fig. [S] left. Therefore we 
have to seek alternative ap proaches. 

In 'Esau ivel et al.l ()2003 ^ it was shown, that if we use 
the thickest possible velocity slice, we get a correct spec- 
trum. Let us consider, why this can happen. Looking 
at the right panel on Fig. [5l we see, that at higher 
wavenumbers ky 's the VCS spectrum is affected by shot- 
noise. As both the VCA and the VCS spectra are just 
different projections of the same 3-d power spectrum in 
the PPV space (see Eqs. [8] and [10]), the noisy har- 
monics seen in the Pi plot are most likely to affect the 
P2 spectrum. If we use only /c^'s, which correspond to 
the intact part of Pi, the picture improves significantly, 
as it can be seen in Fig. [51 left. Such truncation of the 
fct, -domain is equ i valent to taking a thicker velocity slice 
in lEsquivel et al.l (|2003f l. 

It is also important to test the VCS prediction of steep- 
ening of Pi when the beam is made wide enough to emu- 
late the low resolution regime. To do so we used steeper 
velocity spectrum, with slope 4, which is less demanding 
in terms of the number of points along line of sight. If 
we meet condition |15| we have a good agreement with 
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the predictions for both high- and low-resolution regime, 
see Fig. [71 

Besides of testing of analytical results, synthetic data 
calculations could be helpful when pure analytical ap- 
proach is difhcult or impossible. For instance, we 
have tested if VCS technique is applicable in two-phase 
medium such as interstellar HI. It was speculated in 
LPOO that the fluctuations in the channel maps are dom- 
inated by cold (100 K) component of HI, while the warm 
(6000 K) component is subdominant for the VGA stud- 
ies. This corresponds to the decrease of / in Eq. [TU] . 

In our simulations the phases separation was set by a 
threshold, applied to the density field. We assumed that 
the cold phase corresponds to the density larger than 
the threshold. We have also assumed that turbulence is 
the same for both phas es, which can be justified from 
physical arguments fsee iLazarian fc Pogosvanll2000[ ). 

In Fig. [8] the slope of Pi is plotted against the volume 
fraction of emitting medium for different spectral indices 
of velocity. As one can see, with our assumptions there is 
no visible influence of the fraction of the cold phase. We 
argue that the points for the shallow velocity spectrum 
have rather lower quality, because the condition of the 
shot-noise elimination |15| is not met in our simulations. 

5. STUDIES OF ABSORPTION LINES 

We must stress, that although fields of applicability of 
discussed techniques intersect, they do not coincide. It 
is quite evident, that VGA requires less resolution over 
velocity coordinate, and it can be used with lesser Mach 
number than the one needed for VGS. On the other hand, 
even if we have very poor statistics over angular coordi- 
nates and VGA doesn't work, VGS still can be used. It 
allows, for instance, working with absorption lines from 
distant sources, observed through the same turbulent 
medium, like shown on Fig. [9l The limiting case is using 
only one spectral line, but this requires significant reg- 
ular velocity shear (it "unfolds" velocity fluctuations to 
wider spectral line and leads to better statistics, which 
we can utilize by averaging of adjacent /c„-harmonics, see 



Fig. [To]). Without the shear, our simulations show, that 
only 10 independent measurements are needed to gain 
required statistics for N„ 20. 

6. SUMMARY 

We provided numerical testing of the predictions for 
the VGA and the VGS techniques, identified the source 
of shot noise reported in earlier testings, provided an 
analytical estimate of the resolution required to eliminate 
the noise and checked successfully this requirement. 

As the VGA and VGS techniques are closely related 
(cf. Eqs. [8] and [lOj), the presence of shot-noise in VGS 
inevitably affects the quality of VGA spectrum. One has 
to meet a rather demanding requirement given by Eq. 
[TS] to eliminate the noise. 

Shot-noise appears in real observations when the num- 
ber of received photons is low enough. Such could be X- 
ray or optical observations. Gondition |15| gives in this 
case the required number of photons and can be used for 
estimation of observation time. 

We found, that if we accurately deal with shot-noise, 
we get spectra, which are in agreement with predic- 
tions of VGA/ VGS techniques. However the condition 
for shot-noise elimination can be very demanding for 
shallower velocity spectra (Nz — > oo exponentially when 
at, — *■ 3) . This can make the 3-d and even 2-d simulations 
testing for the VGS and the VGA rather difficult. 

We observed, that fragmentation of emitting medium 
due to presence of warm phase does not change the slope 
of Pi, predicted in VGS, which allows using this tech- 
nique for neutral hydrogen in Milky Way in which both 
cold and warm gas are present. 

Having numerical tools to test VGA/VGS techniques 
in various physical environments, we enhanced the relia- 
bility of the results obtained from observational data. 

Acknowledgments Authors thank Dmitry Pogosyan 
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Genter for Magnetic Self-Organization in Laboratory and 
Astrophysical Plasmas. 



APPENDIX 
CODING 

The simulation programs were written in GH — h. There are many reasons for using this language. One of them is 
that we have a convenient infrastructure for managing of input data and initializing of an object hierarchy, provided 
by class data. It gives us an easy way of setting up the input parameters without a necessity to recompile before each 
run. It also allows us to easily configure a batch job and hide details of MPI initialization, if the latter is used. 

The input for the class data is a set of text files with a usual "id=value" syntax, linked to each other in a tree. 
When we call the data constructor with a name of a root file, this tree is mapped to a tree of objects of the type data, 
each representing one input file. 

The object structure of a simulation program is also a tree. Therefore it is convenient to map it to some subtree in 
the data hierarchy. If we do this, constructor signatures rarely become more complex than MyClass : : MyClass (data*) 
and all data initialization is done by instantiating of a root object: 

int main(void) 
{ 

data input ("root . df" ) ; 

DProcessor dp(input .getnode ("DPRDCESSDR") ) ; 
dp.GenerateAndProcessO ; 
dp . SaveResults ; 
return : 



Data files allow references defined by a path with respect to a root node. This allows us, for example, to avoid 
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duplicating of identical data, or to put frequently changed data in one file. Such reorganizations do not require any 
changes in a program because references are transparent for client code. 

This mechanism is also used for running a program in batch mode. In this case one file ( "plugin" ) contains all the 
information to be changed at each run and is re-generated by a controlling program before restarting of simulation. 
A separate node is created for defining this information, a "parameter space" . Practically it means that wc specify 
parameters to be scanned in a parspace.df file and run our runtask command, which takes care of plugin. df, 
controls the execution and stores results of each simulation in a dynamically created directory tree. 

Class diagram for our programs, shown in Fig. [11] is approximately the same regardless of VGA or VCS version, 
number of dimensions, or presence of parallelization. Generation of random field is implemented in an abstract class 
RField, from which classes Velocity and Density are derived. The function RField: : PowerSpectrum(vector k) 
is clean and needs definition in subclasses. This allows us to have different power spectra with the same random 
field generation code. Glass Emissivity implements emissivity law (linear or quadratic) and applies a spatial window 
function. Glass Instrument generates spectral data of required dimension and applies an instrument beam (in 2-d 
and 3-d VGS programs). These data are processed in class PI (VGS) or in class P2 (VGA). Glass DProcessor controls 
generation and processing of realizations and accumulation of spectra. 

The source of simulation programs is available at www. astro. wise. edu/^^lazarian/simulations, see Tab. 2] for 
the complete list. We also provide a utility for running simulations in batch mode, runtask. The programs need 
installed IDL to produce graphics. 
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TABLE 1 
Restoring of velocity spectral 

INDEX av FROM SIMULATIONS 





expected a. 


„ restored a„ 


Nz 




3.56 


3.54±0.01 


8388608 




3.67 


3.65±0.02 


524288 




3.78 


3.79±0.02 


131072 




3.89 


3.89±0.02 


32768 




4.00 


3.91±0.04 






Note. — 


16 realizations 


of spectral 




line have been used for each 








TABLE 2 




VGA PREDICTIONS 


ABOUT P2 SPECTRAL INDEX, 


VELOCITY-DOMINATED MODE 


density spectrum 


2-d picture 


plane spectrum 


1-d picture plane spectrum 


steep 


2 




2 


shallow 


2a5-a„+3 
2 




2 



Note. — See lLazarian fc PoeosvanI li2000l ') for details 



TABLE 3 

VCS PREDICTIONS ABOUT Pi SPECTRAL INDEX, PARALLEL LINES 
OF SIGHT 



density spectrum 


pencil beam 


flat beam 


low resolution 


steep 
shallow 


2 


4 


6 


ap—Z 
2{ae-2) 


c«„-3 
2(ae-l) 


a„-3 
2cie 


c«„-3 


c«„-3 


a„-3 



Note. — See IChepurnov fc LazarianI l|2006l ') for details 



TABLE 4 
List of simulation programs 



program 


dimensions 


techniques 


architecture 


spectld 


1 


VCS 


single processor 


spect2d 


2 


VCS 


single processor 


spect3d 


3 


VCS 


single processor 


p_spect2d 


2 


VCS 


MPI 


p_spect3d 


3 


VCS 


MPI 


vca3d 


3 


VCA 


single processor 


p_vca2d 


2 


VCA 


MPI 


p_vca3d 


3 


VCA 


MPI 


Note. 


— The 


programs 


are available at 



WWW. astro . wise . edu/~lazarian/ simulations 
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PPVCube 




y 



Fig. 1. — Illustration to VGA channel modes. Eddies, exceeding channel width are velocity-dominated, those with extent over v less than 
Av are density-dominated. Velocity-dominated spectrum has the slope {9 — ay) /2 (for steep density spectrum), the density-dominated one 
has the slope equal to (a^ and are the slopes of velocity and density power spectra) 
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Eddie modes: 

1 - low resolution 

2 - high resolution Qioud 

3 - intermediate 




V 



Fig. 2. — Illustration to VCS resolution regimes. Eddies within beam size are in low resolution mode and cause steepening of Pi. The 
ones exceeding the beam size are in high-resolution mode and correspondent Pi is shallower. Therefore high-resolution mode is more 
preferable regarding to the signal-to-noise ratio. 
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Fig. 3. — Illustration to the requirement |15) to the number of points over line of sight. The rows correspond to different velocity spectral 
index av (4, 3.78 and 3.56 from top to bottom). VCS spectra in the middle column correspond to Nz close to the estimation (2^"*, 2^"^ and 
2^^ respectively), in the left column Nz is 4 times lower, in the right — 2 times greater. Undervalued increases shot- noise. (Solid lines 
show the expected slopes). 




1-0 l.T: /,0 3-0 -3-5 -3-0 -2-5 -2-0 -1-5 -1-0 

leg K log kv 



Fig. 4. — Simulations with the number of points over line of sight Nz = 32768, sufficient to suppress shot-noise (Nz = 20000 is needed). 
Left: VGA spatial spectrum (2-d simulation), shallower solid line - expected velocity-dominated spectrum, steeper solid line — density- 
dominated spectrum. Right: correspondent VCS spectrum for high resolution mode (1-d simulation), solid line - expected slope. Velocity 
spectral index is 4. 
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Fig. 5. — Simulations with undervalued number of points over the line of sight {N^ = 32768 instead of 420000). Left: VGA spatial 

spectrum (2-d simulation), shallower solid line - expected velocity-dominated spectrum, steeper solid line - expected density-dominated 
spectrum; the spectrum is distorted (the velocity-dominated one is expected). Right: correspondent VCS spectrum for high resolution 
mode (1-d simulation), solid line - expected slope, shot-noise clearly visible. Velocity spectral index is 11/3, which implies larger minimal 
A^z, than for the case on Fig. |4] 



P2 for dv=375 




Fig. 6. — The same as Fig. O but with noisy fci,-harmonics filtered out when calculating the VGA spectrum (left). 




Fig. 7. — 2-d VCS simulations for high-resolution (left) and low-resolution (right) modes. Solid lines show expected slopes. Velocity 
spectral index is 4. 



Turbulence Spectra from Spectral Lines 



11 





4.5 -| 




- 

4.0- 




- 

3.5- 




3.0- 




2.5- 


o 




(A 


2.0- 






1.5- 




1.0- 




0.5- 




0.0- 



0.01 



high resolution, pencil beam 



-■-a^=3.67 



-I 1 — I — I — I — I — I- 



-I 1 — I — I — h 



0-1 

volume fraction of cool phase 



Fig. 8. — Slope of Pi as a function of an emitting medium volume (for different velocity spectral index a«). 



12 Chepurnov & Lazarian 




Fig. 9. — VCS allows recovery of velocity statistics from absorption lines from stars. The whole Pi is guaranteed to be in high-resolution 
mode, which provides better dynamical range over ky. Numerical simulations show, that very few independent measurements are needed 
to gain required statistics. 
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-'/'^ -?.Ti -■.= -■ ^ -1.? 

One realization, spectaim compression is 16. 
Aif due to galaxy rotation is 8 times greater 
than turbulent line width. 



Fig. 10. — VCS can be applied even to a single spectral line. Configuration can be as shown on this picture. Regular velocity shear from 
the galactic rotation results in statistics being sufficient to recover Pi 



DProcessor 



HGenerateAndProcessC) 
HSaveResultsO 



I 



Instrument 



i-GenerateC) 



I 



P1/P2 



+Accumulate() 

+FinalizeC) 

+SaveResults() 



Emissivity 




Velocity 


+Generate() 

1 


+Generate() 

#PowerSpectrum(k: vector) : double 








Density 




+Generate() 

#PoweTSpectrum(k: vector) : double 









RField 



+GenerateC) 

#PowerSpectrum(k:vector) : double 



Fig. 11. — Class diagram for VCS/VCA simulation programs. Class RField provides functionality for generating of random fields, 
class Instrument provides spectrometric data. Processing class is named PI for VCS and P2 for VCA. Class DProcessor provides coordi- 
nation between Instrument and P1/P2. This diagram is the same for versions with different number of dimensions and with or without 
parallclization. Triangle designates inheritance, rhomb - aggregation. Polymorphic functions arc given in italic. 
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